A Practical Scheme for Error Control using Feedback 
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I. INTRODUCTION 

In the mere space of a decade, quantum information 
theory has blossomed into a burgeoning field of of exper- 
imental and applied research. The initial push for this 
rapid development was provided by Peter Shor's discov- 
ery of an algorithm that enables quantum computers to 
find the period of a periodic function much more effi- 
ciently than any (known) classical computer algorithm 
01 . However, even after Shor's discovery there was much 
doubt about the practicality of quantum computing de- 
vices due to their fragile nature. The coherencies be- 
tween systems carrying the quantum information that are 
crucial to quantum computing algorithms are extremely 
vulnerable and easily destroyed by unavoidable interac- 
tions with the surrounding environment. Furthermore, 
aside from this decoherence, another concern was the ac- 
cumulation of errors introduced by imperfect operations 
performed on the encoded information. 

Both these concerns were largely put to rest by the 
key development of quantum fault tolerance. The error 
accumulation was shown to be tolerable as long as the 
systematic error introduced by each operational element 
was below a critical threshold value jj- This threshold 
result relies heavily upon the development of quantum 
error correction codes. These codes, the first of which 
were discovered by Shor Q and Steane 0|, redundantly 
encode information in a manner that allows one to correct 
errors while preserving coherencies and thus the encoded 
information. 

The main ingredients in the implementation of these 
error control codes are projective von Neumann measure- 
ments that discrctize the errors into a finite set. and fast 



controlled unitary gates that provide the ability to cor- 
rect any corrupted data. Of course, instantaneous pro- 
jective measurements and unitary gates are not perfectly 
implementable in any system; here we will be concerned 
with the details of how one would actually implement er- 
ror correction practically on a system where the physical 
tools available are necessarily physically limited. 

We extend previous work Q and describe an imple- 
mentation of error control that utilizes stabilizer error- 
correcting codes and employs weak measurement and 
Hamiltonian feedback to effectively protect an unknown 
quantum state. This scheme has similarities to the one 
described previously in [3; however, whereas that proto- 
col uses a full state estimation technique that is com- 
putationally intensive, this one uses a simple filtering 
technique that is easily implementable. This protocol is 
therefore a reasonable one for many of today's quantum 
computing architectures. 

This paper is organized as follows. Section IITI reviews 
the key ideas we use: weak measurement, feedback, and 
stabilizer codes. Section IIIII describes our error control 
scheme and two specific instances of it. Section de- 
scribes the simulations performed to analyze the perfor- 
mance of our scheme, presents the results, and considers 
the effect of measurement inefficiency. Section exam- 
ines an actual quantum computer architecture and how 
this error control protocol could be implemented on this 
architecture. Section IVII concludes. 
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II. REVIEW OF CONCEPTS 

A. Quantum feedback control 

In order to describe the behavior of a quantum system 
with feedback, we must first examine the description of 
an open quantum system undergoing continuous weak 
measurement; later we will add feedback conditioned on 
the measurement results. Continuous measurement is 
modelled by considering the system of interest S to be 
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weakly coupled to a reservoir TZ. In order to utilize the 
Born-Markov approximation, we assume that the self- 
correlation time of the reservoir is small compared to 
the time-scales of the system-reservoir coupling and the 
system dynamics. This essentially says that TZ measures 
S continuously but quickly forgets, or dissipates away, 
the result of the measurement. This allows us to write 
the unconditional dynamics of iS as the following master 
equation |^: 

m = -i[H, p{t)] + ^^,.v[cMt)- (1) 

Here p is the reduced density matrix of S, H is the system 
Hamiltonian, {c^} arc the collection of system-reservoir 
interactions (in our case, where we are considering these 
interactions to be weak measurements, these arc the Hcr- 
mitian operators corresponding to the observables) , is 
a parameterization of the coupling strength of , and V 
is the decoherence superoperator given by 

V[A]p = ApA^ ~ i [A'^Ap + pA^A) (2) 

for any operator A. Note that we set h ^ \ throughout 
this paper, except for section Ivl 

This equation describes the unconditional dynamics of 
S because we are assuming that the measurement records 
are ignored. That is, ignoring the measurement records 
corresponding to the observables {c^} means that the 
best description of S we have is one where we average 
over all possible measurement records, and hence all pos- 
sible quantum trajectories that the system could have 
traversed 

If we choose not to ignore the measurement records, 
we instead get a conditional evolution equation for the 
system j^lj: 

m 

dpc{t) = -l[H,p^{t)]dt + K^J)[c^]p^{t)dt 

m 

+ E^^[cM]Pc(t)dW^MW (3) 

where pc is the system density operator conditioned on 
the measurement records of {c^}, dWfi,{t) are Weiner in- 
crements (Gaussian distributed random variables with 
mean zero and autocorrelation {dW{s)dW{t)) = 6{s — 
t)dt) 0], and Ti. is the superoperator 

n[A]p ^ Ap + pA'^ - p ir[Ap + pA^] (4) 

for any operator A. The measurement record from the 
measurement of c^j is: 

dQ^ {t)^Kf, {c^ + cl)cdt + y/l^JWf, (t) (5) 

where (a)^ — tr(pca)- In terms of quantum trajectories, 
this equation corresponds to a diffusive unravelling of the 



master equation given in Eq. From here onwards, 

for simplicity we will specialize to the case of m = 1 in 
equations and (jSJ (i.e. only one c^). 

Given this model, we can now consider adding feedback 
to the system. In general, the feedback will be a function 
of the entire measurement record history. And if we use 
Hamiltonian feedback, the conditional stochastic master 
equation (SME) with feedback becomes 

dpait) = -i[H,pc{t)]dt 

+KD[c]pc{t)dt + y/lm[c]pc{t)dW{t) 
^iRQit)[F, p,it)]dt (6) 

where Rq (t) is some arbitrary function of the entire mea- 
surement history Q{t), and F is the feedback Hamilto- 
nian. Note that all we have done is to add a Hamiltonian 
evolution term whose strength is conditioned by a func- 
tion of the measurement record. 

As shown in 0, in the restricted case that Rqit) is a 
linear function of the measurement value at time t only 
(i.e. RQ{to) = fi inear [Q(io)]), we can simphfy equation 
© further, and derive a master equation for the uncon- 
ditional dynamics of the system. This restricted case 
is often referred to as Markovian feedback, and is con- 
sidered in connection with quantum error control in |^. 
However, in general, when the feedback is conditioned by 
a current that is some arbitrary function of the measure- 
ment history, it is not possible to treat the evolution an- 
alytically, and numerical simulation is the only recourse 
for solving (|SJ). An important special case of this general 
feedback is one in which the function Rgit) is designed 
to compute an estimate of the stat e Pr (t) and output an 
appropriate feedback strength [lol We will refer to 
this as Bayesian feedback, following reference In 
1^ an almost full state estimation is done en route to 
error control, and here we will consider a simpler and 
more practical version of that scheme that performs only 
a partial state estimation. 



B. Stabilizer codes 

In making continuous weak measurements on our sys- 
tem, we would like to choose the measurements in such 
a manner that we gather as much information about the 
errors as possible while disturbing the logical qubits as 
little as possible. These are exactly the conditions satisfi- 
able by encoding the information using a quantum error 
correcting code; and the powerful stabilizer formalism 
|l3l Il4j provides a way to easily characterize many of 
these codes. We will restrict our attention to these stabi- 
lizer codes and in this section provide a brief description 
of the main result of the formalism and give an example. 
For more detailed discussions, the reader is referred to 
[13 and [3. 

We begin by introducing the Pauli group 

p„ = {±i,±z}®{/,x,y,z}®" (7) 
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where X,Y^ and Z denote the Pauh operators (T^,(Ty, 
and (Tj respectively. To simpHfy notation, we wiU omit 
the tensor product symbol when notating members of Pn 
(e.g. ZZI= Z(g)Z(g)I). 

Now, if we encode k logical qubits in n physical qubits 
(n > k), then we can think of errors on our physical sys- 
tem as the action of some subset {Ej} C Pn uM- Thus, 
we would like to choose our encoding in a manner that 
allows us to detect and correct the action of those group 
elements in {Ej}. The main result from the theory of 
stabilizer codes tells us about the possibility of choosing 
such an encoding. It says that provided the elements 
of {Ej} satisfy a certain condition, it is always possi- 
ble to choose a codespace, C, that can be used to detect 
and correct these error elements [T^ . Furthermore, this 
codespace has some special properties: 

1. There exist a set of operators in P„, called the sta- 
bilizer generators and denoted by pi, (72, <7r, such 
that every state in C is an eigenstatc with eigen- 
value -1-1 of all the stabilizer generators. That is, 
gi\ip) = for all i and for all states 1-0) in C. 
Moreover, these stabilizer generators are all mutu- 
ally commuting. 

2. The stabilizer code error correction procedure in- 
volves simultaneously measuring all the stabilizer 
generators and then inferring what correction to 
apply from the measurement results. The formal- 
ism states that the stabilizer measurement results 
indicate a unique correction operation. 

This result tells us that once we identify a set of one 
qubit errors in Pi that we are concerned about, it is pos- 
sible to choose a stabilizer codespace that can be used to 
protect the encoded information against such errors. The 
error detection-correction procedure involves measuring 
each of the stabilizer generators and then applying a cor- 
rection corresponding to the results obtained from the 
stabilizer measurements. 



C. Example: Three qubit bit-flip code 

A common error encountered in quantum computing 
implementations is the bit flip. This type of error has the 
effect of reversing the encoded qubit 's value at random 
times. That is, \'>p) — > Xlip) with probability p, and 
IV') 1^) with probability 1 — p (where again, X = a^, 
and \ip) is a one qubit state). 

One encoding that protects against this type of error 

is 



|0)l 



|000)p 

llll)p 



(8) 



ZZI 


IZZ 


Error 


Correcting unitary 


+1 


+1 


None 


None 


-1 


+1 


on qubit 1 


Xll 


+1 


-1 


on qubit 3 


IIX 


-1 


-1 


on qubit 2 


IXI 



where the right hand side shows the physical encoding in 
three qubits of the logical qubit value on the left hand 
side. That is, C = span{|000), |111)}. The stabilizer 



TABLE I: The three qubit bit-flip code. Note that each error 
results in a diflferent sequence of stabilizer generator measure- 
ment results. 



generators for this codespace are the operators ZZI and 
IZZ. The code can be used to detect and correct any of 
the errors XII, IXI, and IIX. The correction procedure 
involves measuring the two stabilizer generators and then 
applying the appropriate correcting unitary according to 
the rules of table U which corrects for the bit-flip errors. 



III. THE ERROR CORRECTION SCHEME 
A. The general scheme 

Once information is encoded using a quantum error- 
correcting code, conventional error control proposals use 
projective measurements to measure the stabilizer gen- 
erators and fast unitary gates to apply the corrections if 
necessary. In such schemes the detection-correction oper- 
ation, which is initiated by the projective measurement, 
happens at discrete time intervals, and these intervals are 
chosen so that the average number of errors within an 
interval is correctable. We will refer to such implementa- 
tions as discrete error correction schemes because of the 
discrete nature of the detection-correction operation. 

In this section we present a protocol that combines 
weak measurements of the stabilizer generators with feed- 
back to perform continuous error correction. Because of 
the encoding, these measurements will be unobtrusive 
when the system is in the codespace and will give error 
specific information when it is not. However, the require- 
ment of weak measurements makes the measurement cur- 
rents described by © noisy, and therefore ineffective for 
feedback conditioning. In order to use the information 
from these measurements to condition the feedback, we 
must smooth out some of the noise in the currents. The 
smoothing can be easily done using a low-pass filter; how- 
ever, such a filtering process introduces its own compli- 
cations. Specifically, such filtering makes it impossible 
to derive a master equation describing the evolution be- 
cause the noise in the feedback signal at time t is not 
independent of the system state at time t. In essence, 
our smoothing procedure makes Markovian feedback im- 
possible, and leaves the alternative of Baycsian feedback. 

Now, full state estimation is a computationally expen- 
sive procedure - it most often involves solving an SME in 
real-time. In fact, as shown in the resources needed to 
apply a full state estimation feedback procedure to quan- 
tum error control scale exponentially with the number of 
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qubits in the stabihzer code. Fortunately, we do not need 
to do a full state estimation. Instead, the coarse-grained 
state estimate that the stabilizer measurements provide 
— whether the state is in the codespace or not, and if 
not, how to correct back into the codespace — is pre- 
cisely the information needed for error control. That is, 
instead of estimating pdt), we simply need to reliably 
identify the stabilizer generator measurement results (in 
the presence of noise) in order to place Pc{t) inside or 
outside the codespace. Furthermore, as seen in the ex- 
ample of section III CI this information is contained in 
the signatures of the stabilizer generator measurements 
(whether they are plus or minus one), a quantity that is 
fairly robust under the influence of noise. These observa- 
tions suggest that weak measurement and feedback can 
be used to continuously detect and correct errors. 

The general form of the error correcting scheme we 
propose is similar to discrete error control, but with a few 
modifications to deal with the incomplete information 
gained from the weak measurements. The scheme can be 
stated in four steps: 

1. Encode information in a stabilizer code suited to 
the errors of concern. 

2. Continuously perform weak measurements of the 
stabilizer generators, and smooth the measurement 
currents. 

3. Depending on the signatures of the smoothed mea- 
surement currents, form conditioning signals for 
feedback operators on each physical qubit. These 
conditioning currents will be highly non-linear 
functions of the measurement currents because the 
conditional switching based on signatures is a non- 
linear operation. 

4. Apply feedback Hamiltonians to each physical 
qubit, where the strength of the Hamiltonians is 
given by the conditioning signals formed in the pre- 
vious step. 

Given m stabilizer generators and d errors possible on 
our system, the SME describing the evolution of a system 
under this error control scheme is 

d 

k=l 
m 

+ KV[Mi]p,{t)dt + ^U[Mi]p,{t)dWi{t) 
1=1 

d 

+ J2-iGk{t)[Ft,pS)]dt (9) 

k=l 

where jk is the error rate for error Ek, k is the measure- 
ment strength (assumed for simplicity to be the same for 
all measurements Mi), Fk is the feedback Hamiltonian 
correcting for error Ek, and Gk is the feedback condi- 
tioning signal for Fk- Each Gk is a conditional function 



of the signatures of all the smoothed stabilizer measure- 
ments, {Ml}. 

Equation Q has three parts to it: the first line de- 
scribes the effects of the error operators, the second line 
describes the effects of the weak stabilizer generator mea- 
surements, and the third line describes the effect of the 
feedback. [Also note that wc have set the system Hamil- 
tonian, iJ in lO, to zero.] 

This general scheme is illustrated by the following ex- 
amples. The systems described by these examples are 
also the ones simulated in section Hvl 

B. Example: A toy model 

This first example is somewhat artificial, but serves as 
a good illustration of our protocol. The 'codespace' we 
want to protect is simply the state |0), the errors are ran- 
dom applications of X, and the protocol gathers informa- 
tion by measuring the stabilizer generator Z. Obviously 
this 'code' cannot be used for any information process- 
ing, but it is useful for investigating the behaviour of our 
feedback scheme. 

The dynamics of this system before the application of 
feedback are described by the following SME: 

dpc{t) jV[X]pc{t)dt + KV[Z]pc{t)dt 

+yj^H[Z]pcit)dW{t), (10) 

where 7 is the error rate and k is the measurement rate. 
The measurement current has the form 

dQ{t) = 2K{Z)c{t)dt + ^dW{t), (11) 

Now, the measurement of Z reveals whether the sys- 
tems is in the 'codespace' or not because 

Z|0)=+1|0) 

Z|l) = -l|l). (12) 

However, we do not have direct access to {Z)c, but rather 
only to the noisy measurement current Hll|l . Therefore 
we must smooth out the noise on it to obtain error infor- 
mation, and we will choose the following simple filter to 
do so: 

R{t) = ^ ^e-^^'-''^dQ{t') (13) 

This integral is a convolution in time between the mea- 
surement signal and an exponentially decaying signal. 
In frequency space, this acts as a low pass filter, and 
thus the output of this operation is a smoothed version 
of the measurement current with high frequency oscilla- 
tions removed ^ . The filter parameters r and T determine 



This low pass filter is far from ideal. It is possible to design 
low-pass filters with mu ch fi ner frequency selection properties 
(e.g. Butterworth filters) [ifil . and we expect schemes using such 
filters to perform better than this simpler version. 
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the decay rate and length of the filter, respectively, and 
Af = —(1 — e"*"^) serves to normalize R{t) such that it 
is centred around ±1. 

We will use the signature of this smoothed measure- 
ment signal to infer the state of the system and thus to 
condition the feedback. Explicitly, the form of the feed- 
back conditioning current is 

Git)^h^'^ ifi?(t)<o (^4^ 

I otherwise 

Thus, we describe the behaviour of the system with feed- 
back using 

dp^{t) ^ -fV[X]pc{t)dt + KV[Z]pc{t)dt 
+^/^'H[Z]pa{t)dW{t) 
^iXG{t)[X,pc{t)]dt (15) 

where A is the maximum feedback strength. 

Clearly this feedback conditioning current is non- 
Markovian (and non- linear). As mentioned above, this 
makes the Markovian simplification impossible, and 
therefore the most direct route to evaluating this error 
correction protocol is numerical simulation. This is done 
in section Hvl 

We note at this point that there are several open pa- 
rameters in the SME H15|l . These parameters are the 
following: 

1. 7 - This is the error rate and is largely out of the 
experimenter's control. 

2. r - The decay rate of the smoothing filter. Large 
values of r yield responsive measurement cur- 
rents, while small values of r introduce more de- 
lay but make the processed measurement current 
smoother. We expect there to be some optimal 
value of r that achieves a trade-off between respon- 
siveness and smoothing ability. 

r is intimately connected to the other filter pa- 
rameter appearing in T, the size of the fil- 
ter's memory, which determines how many mea- 
surements from the past the filter uses in its cal- 
culations. We will choose T to be large enough 
so that the decaying exponential filter is not trun- 
cated prematurely. A T that is some large enough 
multiple of the filter's time constant, 1/r, would 
be ideal. Since these parameters are dependent on 
each other, we will only consider one of them (r) to 
be free. 

3. A - The maximum strength of the feedback Hamil- 
tonian. The value of this parameter is determined 
by the physical apparatus and the method of feed- 
back. We expect the performance of the protocol 
to improve with A, because increasing A increases 
the range of feedback strengths available. 



4. K - A parametrization of the measurement strength 
used in measuring the stabilizer generators. The 
larger k is the more information we gain from these 
measurements, and thus we expect the performance 
to improve with increasing n. 

In summary, we have three parameters to control - one 
filter parameter, one feedback parameter, and one mea- 
surement parameter. We expect there to be a region in 
this parameter space where this error control scheme will 
perform optimally. We will investigate this issue using 
simulations. 

C. Example: Bit flip correction 

This example is similar to the toy model above but 
looks at a more realistic error control situation. We will 
describe the dynamics of a continuous error correction 
scheme designed to protect against bit flips using the 
three qubit bit flip code of section III CI 

The measurement currents and SME of the system be- 
fore the application of feedback are 

dp,Xt) = -i{V[XII]+V[IXI]+V[IIX])pc{t)dt 
+k{V[ZZI] + V[IZZ])pcit)dt 
+V^iH[ZZI]dWi{t) 

+ n[izz]dW2{t))pc{t) (16) 

dQi{t) = 2K{ZZI)c{t)dt + y/^dWi{t) (17) 
dQ2{t) = 2K{IZZ)c{t)dt + ^dW2it), (18) 

where 7 is the error rate for each qubit, and k is the mea- 
surement strength. We will assume that the errors on 
different qubits are independent and occur at the same 
error rate, and also that the measurement strength is 
the same for both stabilizer generators. (The assump- 
tion of identical rates is made for simplicity and can be 
removed.) 

Now, as detailed in table P the measurements of ZZI 
and IZZ reveal everything about the errors. However, as 
in the toy model, we must smooth the measurement cur- 
rents in order to gain reliable error information. There- 
fore, the steps involved in the error correction scheme are 
the following: 

1. Smooth the measurement currents using the follow- 
ing fllter: 

Rdt) ^Jfj^^ e-''(*"*')dg,(i') j = 1, 2 (19) 

The definition of this filter is analogous to (|13() . 

2. Depending on the signatures of Ri{t) and i?2(i) 
apply the appropriate feedback Hamiltonian. That 
is, 

(a) If Ri{t) < and i?2(i) > 0, apply XII. 

(b) If Ri{t) > and i?2(i) < 0, apply IIX. 
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(c) If Ri{t) < and R2{t) < 0, apply IX I. 

(d) If Ri{t) > and Rait) > 0, do not apply any 
feedback. 

These conditions translate into the following feed- 
back conditioning currents: 

Gi(t) = I ^^^^^ ^^^^^ ° '"^'^ ^'^^^ ^ ° (20) 
I otherwise 

Gsit) = I ^'^^^ ^^^^^ > ° '''''^ ^'^^^ ^ ° (21) 
I otherwise 

G2(t) - I ^'^^^ ° ""''"^ ° (22) 

I otherwise 



turning on feedback. The feedback conditioning signals, 
Gi (t) , can be made non zero as soon as we recognize that 
the smoothed measurement signals arc changing sign. 
That is, the feedback mechanism can be turned on as 
soon as we see a significant shift in the stabilizer gener- 
ator measurements from their error free value: one. We 
can state this 'significant' shift more precisely as a change 
of more than n standard deviations from the mean value 
of one, for n sufficiently large. Thus the choice of n de- 
pends on the signal-to-noise ratio of the smoothed mea- 
surement currents, and therefore on the parameters r and 



IV. SIMULATION RESULTS 



Under this scheme, the SME describing the system dy- 
namics with feedback becomes 

dp^{t) = -/(V[XII]+V[IXI]+V[IIX])pc{t)dt 
+K{V[ZZI\+V[IZZ])pc{t)dt 

+^{n[zzi]dWi{t) + n[izz]dW2{t))pc{t) 

-iX{Gi (<) [XII, p,{t)] + G2 {t) [IX I, p, (t)] 
+ G3{t)[IIX,p,{t)])dt (23) 

where A is the maximum feedback strength, which is as- 
sumed for simplicity to be the same for all the feedback 
Hamiltonians. 

Again, the non-Markovian feedback signals make nu- 
merical simulation the most direct method of solution of 
this SME. 

Also, it is worth noting that even though we gain error 
information from the signatures of the stabilizer gener- 
ators, we do not have to wait until the smoothed mea- 
surement signals, for example ((T^ . fall below zero before 



As a way of evaluating the performance of the general 
error control scheme using weak measurements and feed- 
back, we numerically solved the SMEs described in the 
two examples of section IIIII A comparison of the SMEs 
(list and (|23|l shows that the one qubit toy model has 
all the free parameters of the full three qubit code, and 
therefore is a good model on which to explore the param- 
eter space formed by r. A, and k. This is useful because 
the smaller state space of the toy model makes simulat- 
ing it far more computationally tractable than simulating 
the bit-flip correction example. 



A. The toy model 

We chose to simulate the dynamics of (|15ll by way of 
an associated stochastic Schrodinger equation (SSE) for 
two reasons: (i) it is less computationally intensive, (ii) it 
allows us to look at individual trajectories of the system 
if desired. The form of this associated SSE is as follows: 



d[Mt)) = dN{t){X[Mt))~\Mt))) + V^dW{t){Z-{Z){t)) [Mt)) - -{l-{Z){t)Zf [^,{t))dt - i\G{t)X [Mt))dt 

, m 

I 



where dN{t) is a point process increment (in the number 
of errors) described by 

dN{tf = dN{t) (25) 
E[dN{t)] = jdt. (26) 

That is, dN is a random variable that is either or 1 at 
each time step, and is distributed according to the error 
rate 7. A graph of the process dN{t) would be a sequence 
of Poisson distributed (with parameter 7) spikes. In the 
language of quantum trajectories, this SSE is simply one 
possible unravelling of the SME H15|l . 

The SSE was solved using Euler numerical integration 
with time steps dt = 10~^. When ensemble averages 



were required — that is, when we were interested in the 
behaviour of Pc{t) — 600 trajectories were averaged over. 
To evaluate the performance of the protocol, we used the 
codeword fidelity: F{t) = (?/'(0)|p(t)|V'(0)). Here, |?/'(0)) 
is the initial state of the system, which is taken to be |0) 
unless otherwise specified. 

Figurenshows a sample trajectory from the one qubit 
simulation. The figure shows the expectation value of the 
Z measurement as a function of time and also the super- 
imposed filtered measurement signal, R{t). The transi- 
tions (of expectation value of Z) to —1 are due to errors, 
and the transitions back to +1 are due to feedback cor- 
rection. 
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We used this toy model primarily to gain insight into 
the choice of parameters that lead to optimal error cor- 
rection. The conclusions drawn from exploring the pa- 
rameter space using this one qubit simulation are the 
following: 

1. The decay rate of the filter, r, should be determined 
by the strength of the feedback, A. That is, given 
a strong feedback Hamiltonian, it is necessary to 
have a responsive conditioning current; one with 
little memory. 

2. As expected, the larger the measurement strength 
K, the better the protocol performs. 

3. Performance also improves as A, the feedback 
strength, is increased. This is to be expected be- 
cause as A is increased, the range of the strength 
of the feedback Hamiltonian increases leading to a 
greater degree of control. 

4. The interplay between the two processes - measure- 
ment and feedback - must be considered. In partic- 
ular, larger values of k will yield better performance 
only if these values arc not too much larger than 
the value of A. That is, if the measurement strength 
is much stronger than the feedback strength, the 
measurement process disrupts the feedback correc- 
tion process and makes it ineffective. Therefore the 
magnitude of the measurement strength should be 
less than, or of the same order of magnitude, as the 
feedback Hamiltonian strength. 

Given the strong dependence between parameters that 
these one qubit simulations identify, there are really only 
three free parameters in the system: k, A and 7. Since 
the last is out of the experimenter's control, there remain 
two controllable parameters. In practice, neither of these 
parameters, the measurement strength or the feedback 
strength, are completely configurable. The physical im- 
plementation scheme typically limits the range of these 
parameters, and in section we shall see whether the 
practical ranges for one particular implementation allows 
for error control via this feedback scheme. 

I 



This SSE is of course an unravelhng of SME f^ . and all 
parameters are defined as for that equation. 

As in the one qubit case, we solved this differential 



It is instructive to note that the free parameters of the 
protocol are all physical parameters. That is, the opti- 
mal operating regime of the protocol is defined by the 
system's physical features rather than the introduced fil- 
ter. Therefore, it is possible to design a filter that allows 
the protocol to perform optimally for a given set of phys- 
ical parameters (k and A). 



B. Three qubit code simulation 

Now we move on to the simulation of the three qubit 
bit flip code. This simulation behaves in much the same 
way as the one qubit version, but with one key difference: 
for the one qubit 'code', a double error event - where an 
error occurs on the qubit before we have corrected the 
last error - is not too damaging: in this case, the error 
correcting feedback mechanism detects a traversal back 
into the 'codespace' and thus stops correcting. In the 
three qubit code, this situation is a little more compli- 
cated. Let us consider the situation in which a second 
error happens while a previous error is being corrected. 
If this second error happens to be on the same qubit as 
the one being corrected, then in consonance with the one 
qubit 'code', it is not too damaging. However, if the 
second error is on one of the two qubits not being cor- 
rected, an irrecoverably damaging event occurs, because 
in this case the stabilizer measurements cease to provide 
accurate information about the error location, and the 
protocol's 'corrections' actually introduce errors. 

This identifies a key consideration in any continuous, 
feedback based error correction scheme. The finite dura- 
tion of the detection and correction window means that 
we must choose our parameters with this finite window 
small enough that the probability of an error we cannot 
correct (in this case, two errors on different qubits) is 
negligible. This is analogous to choosing the detection- 
correction intervals in the discrete error control case to 
be small enough to avoid uncorrectable errors. 

The SSE that describes the dynamics of the three qubit 
error correction scheme is 



- mm 



(27) 

I 

equation using an Euler method with timesteps dt = 
10~^. Again, ensemble averages were done over 600 tra- 
jectories when needed. The initial state used was |000), 



d\Mt)) = dNi{t){XII\Mt))-\Mt))) + dN2{t){IXI\Mt))~\Mt))) + dN3{t){IIX\Mt)) 

+ y/^ dWi{t){ZZI - {ZZI){t)) \lPa{t)) + y/^ dW2{t){IZZ ^ {IZZ){t)) l^cit)) 

- ^il-{ZZI)it)ZZI)^ \Mt))dt - ^{l-{IZZ){t)IZZf \Mt))dt 

- iXGi(t)XII \4>c{t))dt - iXG2{t)IXI \Mt))dt - i\Gi{t)IIX \4>c{t))dt. 
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and the performance was measured using the codeword 
fidehty F3(^) = (000|p(t)|000). A true Mehty measure of 
the protocol performance would average over all possible 
input states; however, because |000) is most susceptible 
to bit flip errors, the fldelity we use can be considered a 
worst case performance analysis. 

The performance of the error correction scheme using 
this code is summarized by figure |21 This figure shows 
the fidelity versus time curves (F^it)) for several values 
of error rate (7). Each plot also shows the fidelity curve 
(Fi (t)) for one qubit in the absence of error correction. A 
comparison of these two curves shows that the fidelity is 
preserved for a longer period of time by the error correc- 
tion scheme for small enough error rates. Furthermore, 
for small error rates (7 < 0.3) the ^3(i) curve shows a 
vast imrpovcment over the exponential decay in the ab- 
sence of error correction. However, we see that past a 
certain threshold error rate, the fidelity decay even in 
the presence of error correction behaves exponentially, 
and the two curves look very similar; past the threshold, 
the error correcting scheme becomes unable to handle 
the errors and becomes ineffective. In fact, although not 
completely evident from the figure, well above the thresh- 
old the performance of the scheme becomes worse than 
the unprotected qubit's performance. This poor perfor- 
mance results from the feedback 'corrections' being so 
inaccurate that the feedback mechanism effectively in- 
creases the error rate. 

The third line in the plots of figure |21 is of the average 
fidelity achievable by discrete quantum error correction - 
using the same three qubit code- when the time between 
the detection-correction operations is t. The value of this 
fidelity {F^dit)) as a function of time was analytically 
calculated in 



F3, = i(2 + 3e-^^*-e-6^*). 



(28) 



A comparison between i^3(i) and F^dit) highlights the 
relative merits of the two schemes. The fact that the 
two curves cross each other for large t indicates that 
if the time between applications of discrete error cor- 
rection is sufficiently large, then a continuous protocol 
will preserve fidelity better than a corresponding discrete 
scheme. In fact, this comparison suggests that a hybrid 
scheme, where discrete error correction is performed rel- 
atively infrequently on a system continuously protected 
by a feedback protocol, might be a viable approach to 
error control. 

All the F3 (t) curves show an exponential decay at very 
early times, t w to t w 0.1. This is an artifact of the 
finite filter length and our specific implementation of the 
protocol. In particular, our simulation does not smooth 
the measurement signal until enough time has passed to 
get a full buffer of measurements; that is, filtering and 
feedback only start a.t t = T. Of course, this can be 
remedied by a more complicated scheme that smoothes 
the measurement signal and applies feedback even when 
it has access to fewer than T/dt measurements. 
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FIG. 1: A sample trajectory of the one qubit "code" with 
feedback. The top graph just shows the expectation value of 
Z, and the bottom graph shows expectation value of Z and 
the filtered signal R{t). 
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FIG. 2: Fidelity curves with and without error correction for 
several error rates. The thick solid curve is the fidelity of 
the three qubit code with error correction, F3{t) (parameters 
used: dt = 10"*s, k = 150 Hz,\ = 150 Hz,r = 20 Hz,T = 
1500 X dt s). The dotted curve is the fidelity of one qubit 
without error correction, Fi{t). And the thin solid curve is 
the fidelity achievable by discrete quantum error correction 
when the duration between applications is t, F^dit)- 



C. Inefficient measurement 

We have modeled all our measurement processes as 
being perfect. In reality, detectors will be inefficient and 
thus yield imperfect measurement results. This ineffiency 
is typically represented by a parameter rj that can range 
from to 1, where 1 denotes a perfect detector. How 
is this feedback protocol affected by non-unit efficiency 
detection? 
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To examine this question, we simulated the three qubit 
code with inefficient detection. Tlie evolution SME and 
the measurement currents in the presence of inefficient 
detection are as follows: 

dpc{t) = -/{V[XII]+V[IXI]+V[IIX])pc{t)dt 
+k{V[ZZI] + V[IZZ])pcit)dt 
+y/l^{H[ZZI]dWi{t) + H[IZZ]dW2{t))pc{t) 
-i\{Gi{t)[XII, p,{t)] + G2{t)[IXI, p,it)] 
+ G^{t)[IIX,p,{t)])dt (29) 
dQi{t) = 2K^{ZZI)c{t)dt + y/^dWi{t) (30) 
dQ2{t) = 2n^{IZZ)^{t)dt + ^dW2{t) (31) 

where < 77 < 1 is the measurement efhciency, and all 
other quantities arc the same as in equations H16|l and 

The results of these simulations are summarized by 
figure 121 The decay of fidelity with decreasing 77 indi- 
cates that inefficient measurements have a negative effect 
on the performance of the protocol as expected. How- 
ever, the slope of the decay is very small — in particular, 
the graph does not exponentially decay as do Markovian 
feedback protocols — and this suggests that this protocol 
has a certain tolerance to inefficiencies in measurement. 
This is reasonable because the filtering process in the 
protocol has the effect of improving the quality of the 
measurements and thus negating some of the ill effects 
of the inefficient measurements. Also, as in the full state 
estimation protocol of 0, because the feedback condi- 
tioning current is a function of a measurement record 
history — as opposed to just the current measurement — 
errors induced by inefficient measurement tend not to be 
so damaging. Here we see the true strength of this error 
correction scheme: it combines the robustness of a state 
estimation based feedback protocol with the practicality 
of a Markovian feedback protocol. 

V. LINKS TO EXPERIMENT 

In this section we study the possibility of applying this 
error correction technique to a particular quantum com- 
puting architecture. 

A. Solid-state quantum computing with RF-SET 
readout 

Several schemes for solid-state quantum computing 
have been proposed These use the charge 

or spin degree of freedom of single particles to represent 
logical qubits, and measurement involves probing this de- 
gree of freedom. 

Here we examine the weak measurement of one such 
proposal that uses coherently coupled quantum dots 
(CQDsj and an electron that tunnels between the 
dots|23|. The dots arc formed by two P donors in Si, 
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FIG. 3: Average fidelity after a fixed amount of time as a 
function of 1-efficiency for several error rates (parameters 
used: dt = 10~*s, k = 50 Hz,\ = 50 Hz.r = 10 Hz,T = 
1500 xdt s). 

separated by a distance of about 50nm. Surface gates are 
used to remove one electron from the double donor sys- 
tem leaving a single electron on the P-P+ system. This 
system can be regarded as a double well potential. Sur- 
face gates can then be used to control the barrier between 
the wells as well as the relative depth of the two wells. 
Using surface gates, the wells can be biased so that the 
electron can be well localised on either the left \L) or the 
right of the barrier. These (almost) orthogonal lo- 
calised states are taken as the logical basis for the qubit, 
|0) = |1) = \R). It is possible to design the double 

well system so that, when the well depths are equal, there 
arc only two energy eigenstatcs below the barrier. These 
states are the symmetric ground state |-|-) and the an- 
tisymmetric first excited state |— ). A state localised on 
the left (right) of the barrier is then well approximated 
as a linear superposition of these two states, 

\L) = -^(|+) + |-)) (32) 

\R) = ^(1+) - |-» (33) 

An initial state localised in one well will then tunnel to 
the other well at the frequency A = {e^ — e-)/fi where 
e± are the two energy eigenstatcs below the barrier. 

The Pauli matrix, Z = \L){L\ — \R){R\, is diagonal in 
this localised state basis. The Hamiltonian for the system 
can be well approximated by 

H = H^Z + H^X (34) 

where X = \L){R\ + \R){L\. Surface gates control the 
relative well depth huj{t) (a bias gate control) and the 
tunnelling rate A(t), (a barrier gate control) which are 
therefore time dependent. For non-zero bias the energy 
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gap betwee n the ground s tate and the first excited state is 
E(t) = h^Lu{tY + ^{iY- Further details on the validity 
of this Hamiltonian and how well it can be realised in the 
PP"'' in Si system can be found in [2l] |. 

A number of authors have discussed the sources of de- 
coherence in a charge qubit system such as this onej2tl 
I21I . For appropriate donor separation, phonons can 
be neglected as a source of decoherence. The dominant 
sources of decoherence then arise from fluctuations in 
voltages on the surface gates controlling the Hamilto- 
nian and electrons moving in and out of trap states in 
the vicinity of the dot. This latter source of decoherence 
is expected to occur on a longer time scale and is largely 
responsible for 1// noise in these systems. In any case 
both sources of decoherence can be modelled using the 
well known spin-boson model The key element of 

this model for the discussion here is that the interaction 
energy between the qubit and the reservoir is a function 
of ^. 

If the tunnelling term proportional to A(i)X in Eq. 
(|34|l were not present, decoherence of this kind would 
lead to pure dephasing. However, in a general single 
qubit gate operation, both dephasing and bit flip errors 
can arise in the spin-boson model. We can thus use the 
decoherence rate calculated for this model as the bit flip 
error rate in our feedback error correction model. We will 
use the result from the detailed model of HoUenberg et 
al. for a device operating at lOK, and set the error 
rate 7 = 1.4 x 10^s~^. This rate could be made a factor 
of ten smaller by operating at lower temperatures and 
improving the electronics controlling the gates. 

We now turn to estimating the measurement rate , k, 
for the PP+ system. In order to readout the qubit in the 
logical basis we need to distinguish a single electron in the 
left or the right well quickly and with high probability of 
success (efficiency). The technique of choice is currently 
based on radio frequency single electron transistors (RF- 
SET)|23|- We will use the twin SET implementation of 
Buehler et al. IH. 

In an RF SET the Ohmic load in a tuned tank circuit 
comprises a single electron transistor with the qubit act- 
ing as a gate bias. The two different charge states of the 
qubit provide two different bias conditions for the SET, 
producing two different resistive loads, and thus two lev- 
els of power transmitted through the tank circuit. The 
electronic signal carries a number of noise components: 
for example, the Johnson-Nyquist noise of the circuit, 
random changes in the SET bias conditions due to fluc- 
tuating trap states in the SET, etc. The measurement 
must be operated in such a way that the charge state 
of the qubit can be quickly discerned as a departure of 
the signal from some fiducial setting, despite the noise. 
Clearly it takes some minimum time interval, Im, to dis- 
criminate a qubit signal change from a random noisy 
fluctuation. We need to keep the measurement time as 
short as possible. However if the measurement time is 
too short, one may mistake a large fluctuation, due to a 
non-qubit based change in bias conditions, for the real 



signal. In other words one may mistake a 1 for a 0, 
and vice versa. The probability of this happening is the 
efficiency of the measurement, 77 (i^/), which depends on 
the measurement time. The key performance parameters 
are (i) the measurement time, tM^ and (ii) the efficiency 
rjitj^i). An additional parameter that is often quoted is 
the minimum charge sensitivity per root hertz, S. Given 
S determines a minimum change in the charge, A(7, 
that can be seen by the RF-SET at a given bias condi- 
tion. In a measurement time of Im = 6 x lO^^s 
was found for a signal of Aq = 0.2e and an efficiency of 
10~^. We now need to relate this measurement time to 
the measurement decoherence rate parameter, k, of our 
ideal feedback model. 

If the measurement were truly quantum limited (that is 
to say, the signal-to-noise ratio is determined only by the 
decoherence rate k) , the inverse measurement time would 
be of the same order of magnitude as the decoherence 
rate (see |2^). The measurement described in Buehler 
et al.|25j| will almost certainly not be quantum limited. 
However, here we will assume the measurement to be 
quantum limited, so as to obtain a lower limit to the mea- 
surement decoherence rate. Thus we take n = lO^s^^. 

We next need to estimate typical values for the feed- 
back strength. From Eq. H15|l we see that the feedback 
Hamiltonian is proportional to an X operator. In the 
charge qubit example, this corresponds to changing the 
tunneling rate for each of the double dot systems that 
comprise each qubit. The biggest tunneling rate (A) oc- 
curs when the bias of the double wells makes it symmet- 
ric. In |2ll |. the maximum tunneling rate is about 10^ 
s~^, for a donor separation of 40nm. A large tunneling 
rate makes for a fast gate, and thus a fast correction op- 
eration. Thus the maximum value of A can be taken to 
be 10^ s-i. 

To summarise, in the PP+ based charge qubit, with 
RF-SET readout, we have 7 « k « lO^s"^, and A « 
IQSs-i. 

The fact that the measurement strength and the er- 
ror rate are of the same order of magnitude for this ar- 
chitecture is a problem for our error correction scheme. 
This means that the rate at which we gain information is 
about the same as the rate at which errors happen, and 
it is difficult to operate a feedback correction protocol 
in such a regime. Although it is unlikely that the mea- 
surement rate could be made significantly larger in the 
near future, as mentioned above it is possible that the 
error rate could be made smaller by improvements in the 
controlling electronics. Thus it is interesting to consider 
how low the error rate would have to be pushed before 
our error control scheme becomes effective. To answer 
this question we ran the three qubit bit flip code sim- 
ulation using the parameters stated above and lowered 
the error rate until the error control performance was ac- 
ceptable. We found that the fidelity after 1ms could be 
kept above 0.8 on average if the error rate, 7, is below 
10^ s-i (with K = 10^ s-i, and A = 10'^ s"!). So we see 
that a difference in order of magnitude of four between 
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the measurement and feedback strengths, and the error 
rate, is about what this protocol (using the three qubit 
code) requires for reasonable performance. That is, we 
require 

- « - « 10^ (35) 

7 7 

Of course, depending on the performance requirements 
this ratio may be larger or smaller. Also, a full optimiza- 
tion of the filter used in the scheme is likely to drive this 
ratio down by up to an order of magnitude. 

We can compare the requirements of the three-qubit 
code with the one-qubit version. Given the same mea- 
surement and feedback parameters {k = 10^ s~^, A = 10^ 
s~^), the one-qubit 'code' can keep the fidelity above 0.8 
after 1ms when K/7 « A/7 « 10. That is, only one or- 
der of magnitude difference is required between the error 
rate and the measurement and feedback rates. This sug- 
gests that a key issue with feedback based error correction 
schemes is scalability. The ratio between measurement 
and feedback rates and error rate has to increase along 
with the error correcting code size (in qubits). 

VI. DISCUSSION AND CONCLUSION 

We have described a practical scheme for implement- 
ing error correction using continuous measurement and 
Hamiltonian feedback. We have demonstrated the va- 
lidity of the scheme by simulating it for a simple error 
correction scenario. 

As the simulations show, this error control scheme 
can be made very effective if the operational parameters 
(measurement strength, feedback strength, filter param- 
eters) are well matched to the error rate of a given sytem. 
At the same time, the scheme uses relatively modest re- 
sources and thus is easy to implement, as well as being 
robust in the face of measurement inefficiencies. 

From a quantum control perspective, an interesting 
feature of this protocol is the encoding. That is, despite 
using state estimate feedback, the protocol requires little 
side processing due to the fact that instead of a full state 
estimate, it uses a coarse-grained state estimate natu- 
rally suggested by the encoding. In control theory terms, 
this simplification is a result of the specific choice of con- 
trol state space (what to observe and control); a choice 



dictated by the stabilizer encoding and measurements. 
It would be interesting to examine the general conditions 
under which an encoding is available that allows for prac- 
tical, efficient state estimate feedback control. 

The possibility of using continuous error correction in 
combination with its discrete counterpart is an interest- 
ing possibility. Such a scheme has the potential to signif- 
icantly improve the stability of quantum memories, and 
the implications of such a combination scheme for fault 
tolerance would be worth investigating. 

We also studied a solid-state quantum computing ar- 
chitecture with RF-SET readout and the feasibility of 
implementing this error correction protocol on it. Al- 
though the measurement and feedback rates currently 
possible on this architecture do not allow for error cor- 
rection via this feedback scheme with the intrinsic error 
rate, it is foreseeable that as the controlling technology 
improves, this error control scheme will become possi- 
ble on this architecture. From numerical simulations, we 
found the approximate parameter regime where the three 
qubit code using this scheme becomes effective - that is, 
exactly how much improvement is necessary before the 
scheme becomes feasible. It would be interesting to in- 
vestigate this further and explore more rigorously how 
values of k/7 and A/7 dictate protocol performance as 
well as the exact dependency of these parameter ratios 
on the code size. Such an investigation will be crucial 
in addressing the issue of scalability of this error control 
scheme. 
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